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ABSTRACT 


r 

We  say  tha^t  the  random  variable  X  is  more  variable  than  Y  if 
E[f(X)]  >^E[f(Y)]  for  all  increasing  convex  functions  f  .  We  prove 
a  preservation,  under  random  sized  sums,  property  of  this  ordering 
and  then  apply  it  to  branching  processes  and  shock  models. 


SOME  APPLICATIONS  OF  A  RESULT  CONCERNING  VARIABILITY  ORDERINGS 

by 

Sheldon  M.  Ross 


1.  A  VARIABILITY  RESULT 

If  and  X2  are  random  variables  having  respective  distributions 

and  F^  ,  then  we  say  that  X^^  £  X2  (read  X^^  is  less  variable  than  X2) 

or  equivalently  that  F  <_  F  if 

V  ^ 


00  00 
J  f(x)dF^(x)  1 J  f(x)dF2(x) 


for  all  increasing  convex  functions  f  .  Some  easily  derived  properties  of 
this  ordering  are 


1.  F  £  F-  if  and  only  if 

i  Y  ^ 


00  00 

J  Fj^(x)dx  i 


(x)dx  £  I  F2(x)dx  for  all  a 


where  F.  =  1  -  F.  . 
i  1 


2.  If  F  £  G  ,  i  =  1  ,  2  then  F  *  F„  £  G  *  G„  where  *  denotes 

XyX 


convolution. 


We  will  now  present  a  theorem  concerning  this  orderi.'.g  and  in  Sections 


2  and  3  apply  it  to  branching  processes  and  shock  models. 


Theorem  1: 


Let  Xj^,X2,  ...  be  a  sequence  of  nonnegative  independent  and  identically 
distributed  random  variables  and  similarly  '''  '  ^ 

integer  valued  nonnegative  random  variables  that  are  independent  of  the  X^ 
and  sequences.  Then 

N  M 

X  >  Y  i  >  1  ,  N  >  M  -o  I  X  >  I  Y  . 

V  ^  V  i=l  ^  i=l  ^ 


Proof ; 

We  will  first  show  that 


N 

I 


i“l 


M 

Xi>  I 

V  i=l 


Let  h  denote  an  increasing  convex  function.  To  prove  the  above  we  must 
show  that 


'  1^  \ 

(1) 

E 

il 

Since  N  ^  M  ,  and  they  are  independent  of  the  X.  ,  the  above  will  follow 
v  ^ 

if  we  can  show  that  the  function  g(n)  ,  defined  by 

g(n)  -  E[h(X,  +  ...  +  X  )] 
i  n 

is  an  increasing  convex  function  of  n  .  As  it  is  clearly  increasing 
since  h  is  and  each  X^  is  nonnegative  it  remains  to  show  that  g  is 
convex,  or,  equivalently,  that 


(2) 


g(n  +  l)  -  g(n)  is  increasing  in  n  . 
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To  prove  this  let  ^  ^1  ’  i^ote  that 


g(n  +  1)  -  g(n)  -  E[h(S^  +  -  h(S^)] 


E[h(S^  +  -  h(S^)  I  -  t]  -  E[h(t  +  -  h(t)] 

»  f(t)  (say). 


As  h  is  convex,  it  follows  that  f(t)  is  increasing  in  t  .  Also,  as 

S  increases  in  n  ,  we  see  that  E[f(S  )]  increases  in  n  .  But 

n  n 


E[f(S^)]  “  g(n  +  1)  -  g(n) 


and  thus  (2)  and  (1)  are  satisfied. 
We  have  thus  proven  that 


N  M 

I  X  1  I  X 

1  V  1  ^ 


and  the  proof  will  be  completed  by  showing  that 


M  M 

I  X  1  I  Y 

1  V  1 


or,  equivalently,  that  for  increasing,  convex  h 
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2.  A  BRANCHING  PROCESS  APPLICATION 

Consider  two  Galton  Matson  branching  processes  in  which  individuals 

at  the  end  of  their  lifetime  give  birth  to  a  random  number  of  offspring. 

Let  j  >  1  ,  n  >  0  denote  the  number  of  offspring  of  the 

jn  —  - 

th  th 

individual  of  the  n  generation  in  the  i  branching  process,  i  =  1  ,  2  . 

Suppose  that  the  random  variables  ,  j  >  1  ,  n  >  0  are  independent 

jn  -  - 

for  i  =  1  ,  2  and  have  a  distribution  not  depending  on  j  .  In  addition, 
suppose  that 


x(l)  ^  for  all  n  ,  j  . 

-1  Tl  —  1  «  »  J 


V 


/ « \  til  t h 

Let  Z'  ,  i  =  1  ,  2  denote  the  size  of  the  n  generation  of  the  i 
n 

process. 


Proposition  2: 
If  = 


1  ,  i  =  1  ,  2  ,  then  for  all  n 

n  —  n 

V 


Proof; 

The  proof  is  by  induction  on  n  .  As  it  is  true  for  n  =  0  , 
assume  it  for  n  .  Now 


z<i> 

z^i;  =  !  x^^> 

n+l  3,n 


(2) 

n+l 


z<2) 

n 

I 

j-i 


and  so  the  result  follows  from  Theorem  1. | ] 
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We  now  show  that  if  the  second  (less)  variable  process  has  the  same 
mean  number  of  offspring  per  individual  as  does  the  first  then  it  is 
less  likely,  at  each  generation  to  become  extinct. 


Corollary  3 : 


Suppose 

■  i 

X 

1 _ 1 

and  >  y. 

for  all 

j  .  ' 

in  — 

■'  V 

jn 

p{z<l>  .  o} 

If  =  1  .  i  =  1  ,  2 


Proof : 


From  Proposition  2  we  have  that  >  Z^^^  and  thus 

n  -  n 


or,  equivalently,  since 


i  4  i  p  -  1 1} 


which  proves  Che  result. |1 
Remarks : 

(i)  In  [2],  Freedman  and  Purves  showed  that  among  all  branching 

processes  for  which  P[X.  =  1]  =  0  and  E[X.  ]  =  M  <  2  , 

all  j  ,  n  ,  the  one  having  the  least  chance  of  going  extinct 

is  the  one  with  P{X.  =  0}  =  1  -  M/2  =  1  -  P{X.  =  2}  . 

in  jn 

This  also  follows  from  Corollary  3  upon  application  of  the 
following  lemma  (with  o  •  0)  . 
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Lemma  4: 

Let  PCX  =  1}  =  a  ,  PCX  =  0}  =  (1  -  a)  -  ,  PCX 

and  let  Y  be  a  nonnegative,  integer  valued,  and  such  that 
and  E[Y]  =  M.  If  a<M<2-a  then  X  Y  . 

V 

Proof : 

We  must  show  that 

00  00 

I  PCX  1  i}  I  PCY  i}  ,  n  =  1,2,  . . .  . 
i=n+l  i=n+l 

As  E[X]  =  E[Y]  =  M  this  is  equivalent  to 

n  n 

I  ^  i}  i  I  P{Y  >  i)  ,  n  =  1,2,  ...  . 
i=l  i=l 

When  n  *  1  the  above  reduces  to  PCX  =  0}  <_ PCY  *  0}.  This  follows 
since,  as  PCY  =  1}  1  PCX  =  1}  ,  if  PCY  =  0}  <  PCX  =  0}  then  it  would 
not  be  possible  for  E[Y]  to  equal  E[X]  .  When  n  >  1  the  above  is 
equivalent  to 

n 

M  ^  I  PCY  ^  i} 
i=l 

which  follows  since  E[Y]  =  M  . | [ 


PCY  =  1}  <  a 
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3.  A  SHOCK  MODEL  APPLICATION 

Suppose  that  shocks  occur  in  accordance  with  a  renewal  process  having 
interarrival  distribution  G  and  mean  .  Each  shock  gives  rise  to  a 

nonnegative  random  damage  which,  independent  of  all  else,  has  probability 
distribution  F  .  The  damages  are  assumed  to  be  additive  and  we  let 
D(t)  denote  the  damage  at  time  t  .  That  is 


N(C) 


D(t)  =  1  X. 

i=l  ^ 


where  is  the  damage  of  the  i^*^  shock  and  N(t)  is  the  number  of 

shocks  by  t  .  The  system  is  assumed  to  fail  the  first  time  that  D(t) 
exceeds  some  constant  c  .  That  is,  the  system  fails  at  time  T_  „  where 


T_  ^  =  min  { t  :  D(t)  >  c} 

r  ,ti 


We  will  obtain  a  variabilitv  result  about  T  when  both  F  and  G 

F  ,G 

are  NBUE  distributions,  where  a  distribution  of  a  nonnegative  random 
variable  X  is  said  to  be  NBUE  (new  better  than  used  in  expectation)  if 

E[X  -  t  I  X  >_  t]  ^  E(X]  for  all  t  ^  0  . 

Letting 


N(c)  =  max  {n:X,  +...+X  <c} 

1  n  — 


then  the  system  will  fail  at  the  time  of  the  N(c)  +  1  shock. 


Lemma  5: 


If  F  is  NBUE  then 


N(c)  +  1  <  N  (c)  +  1 


v 
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where  N  (c)  is  a  Poisson  random  variable  with  mean  ciJ_  where 

F 

=  E[X]  . 


Proof : 

As  N(c)  is  just  the  number  of  renewals  by  time  c  of  a  renewal 
process  whose  interarrival  distribution  is  NBUE  with  mean  u  the  result 
follows  from  Theorem  3.17  on  page  173  of  [1].|| 

Proposition  6: 

If  F  and  G  are  both  NBUE  distributions  then 


T  <  T 


where  E^  and  E^  are  exponential  random  variables  having  the  same  means 
as  F  and  G  respectively. 

Proof; 

We  can  express  T„  „  by 


T 


F.G 


N(c)+1 

I 

i=l 


Y. 

1 


where  the  Y^  ,  i  ^  1  ,  are  the  interarrival  times  between  successive 
shocks.  They  are  thus  independent  and  have  distribution  G  .  Now  it  is 
well  known  that  an  NBUE  distribution  G  is  less  variable  than  an  exponential 
distribution  with  the  same  mean  and  so 

Y_,  <  e,  when  e.  is  exponential  with  mean  vi„  . 
i  —  i  1  u 


The  result  now  follows  from  Lemma  5  and  Theorem  1. 
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Remark: 


N*(c)+1 


As  T  „  =  T  e  ,  it  follows  upon  conditioning  on  N  (c)  that 

^l’^2  i=l  ^ 


i  )  «  -PC  (UpC)^ 


Gi^l(x) 


where  C  (x)  is  the  gamma  distribution  with  parameters  n  and  l/u^ 
n  ij 

(its  mean  is  np  )  .  Also  if  F  and  G  are  NBUE  then  from  Proposition  6 

(j 

all  of  the  moments  of  T  are  no  greater  than  the  corresponding  moments 

t  ,& 

of  T_  „  .  For  instance 

^l’^2 


"  E[(N*(c)  +  l)Pg]  =  (cPy  +  1)1 


G  • 
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